O 
O 

(N 



Appearance of quark-hadron duality in the Rein-Sehgal model 

Krzysztof M. Graczyk, Cezary Juszczak, Jan T. Sobczyk a * 

a Institute of Theoretical Physics 

University of Wroclaw 

pi. M. Borna 9, 50-204, Wroclaw, Poland 



(N 
> 

i> 
o 

o 

O 

Oh! 



X 



Quark-hadron duality in neutrino-nucleon reactions is investigated under the assumption that cross sections 
in the resonance region are given by the Rein-Sehgal model. The quantitative analysis of the duality is done 
by means of appropriate integrals of the structure functions in the Nachtmann variable. We conclude that with 
the definition of the resonance region W G (M + GeV) the duality holds for neutrino-proton reaction F2 

structure function for Q 2 G (0.5, 3) GeV 2 and it is absent for neutrino-neutron reaction. 



1. INTRODUCTION 

The discovery of quark-hadron duality in 
electron-nucleon interactions £Q rises a natural 
question if the same phenomenon can be seen also 
in neutrino-nucleon scattering. Unfortunately the 
available data is not yet precise enough to dis- 
cuss the problem on the experimental level. At 
present the only possibility is to analyze the ex- 
isting theoretical models of resonance production. 
This approach was adopted by Sato and Lee |2j 
who discovered that for their model of A produc- 
tion the quark-hadron duality is seen: the reso- 
nance peaks of structure functions calculated at 
Q 2 = 0.4, 1, 2, 4 GeV 2 slide along the DIS struc- 
ture functions (with CTEQ6 PDFs) calculated at 
10 GeV 2 , both as functions of the Nachtmann 
variable. 

In this contribution we wish to report an in- 
vestigation of the Rein-Sehgal (RS) model [3] of 
resonance production. The RS model is used by 
almost all Monte Carlo generators of neutrino 
interactions and understanding its properties is 
of practical value. The model describes reso- 
nance production in the region of hadronic in- 
variant mass W below 2 GeV by summing con- 
tributions from 18 resonances. It contains also 
a non-resonant background fine tuned in order to 
get a good agreement with the existing single pion 
production (SPP) data. 



We calculate the structure functions Fi,F2, F3 
as they are defined by the RS model 0]. They 
turn out to be linear combinations of cross sec- 
tions for three polarization states of the inter- 
mediate W boson. The aim of the RS model is 
to describe only SPP channels and the elasticity 
coefficients are used in order to select exclusive 
SPP channels from the overall resonance produc- 
tion cross section. The non-resonant background 
is adjusted only for SPP channels. The analysis 
of duality requires the knowledge of the inclusive 
cross section in the resonance region. In the spirit 
of the RS model one should add contributions 
from inelastic channels taking care of complicated 
interference patterns. In order to make the anal- 
ysis simpler we introduce lir functions which are 
defined as probabilities that in the given region 
of the kinematically allowed space the final state 
is that of SPP. We calculate them numerically 
for each SPP channel separately using the LUND 
fragmentation and hadronization routines. All 
the details of our approach can be found in @j. 

2. FORMALISM AND RESULTS 

The lir functions are defined (for each exclusive 
channel separately) as follows : 
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Figure 1. A comparison of lir functions defined in 
Eq. with resonance elasticity factors taken from 
PDG 5]. In the case of v- neutron scattering we show 
the sum of two functions corresponding to v^n — > 
fj,~pir° and v^n — > (i~rnr + reactions 



It turns out that fi v (W, Q 2 ) do not depend on 
Q 2 . Their plots for CC vN channels are pre- 
sented in Fig. ^ In the case of neutron we show 
the sum of l7r functions corresponding to two ex- 
clusive SPP channels. On the same plot we show 
also the available data on elasticity of N and A 
resonances. We see that Itt functions provide sat- 
isfactory average description. 

We re-scale RS model structure functions and 
obtain approximations of structure functions in 
the RES region: 



l!ES (x Q 2 ) = F t RS ( x >Q 2 *) 

fl7r(x,Q 2 ) 



i = 1,2,3. 



(2) 



In our quantitative analysis we define the fol- 
lowing function: 



n 2 {Q 2 -Q 2 DIS ) 



dxF 2 {Z,Q 2 



(3) 



d£F 2 DIS (Z,Ql IS ) 



£ is the Nachtmann variable. (, m ax.min ar e de- 
fined by the conditions W{^ max , Q 2 ) = M + 
W(£ m in,Q 2 ) = W max . We investigate three 



Figure 2. Double differential cross section ^dQ' 1 f° r 
neutrino energy 3 GeV. 



choices for W max . The perfect duality would 
mean that U 2 (Q 2 ;Q 2 DIS ) — 1- 

We select the physically relevant region in the 
Q 2 variable as Q 2 < 3 GeV 2 by looking at double 

.2 

differential cross section d ^j^ 2 for neutrino en- 
ergy in a few GeV region, see Fig. 2 6 . Having in 
mind MC generators of events the quark-hadron 
duality should ensure the smooth passage of the 
cross section for RES and DIS contributions in 
the relevant kinematical transition region. 

In Fig.|21we show F 2 structure functions for CC 
reactions defined by the RS model for three values 
of Q 2 : 0.4, 1, 2 GeV 2 and also DIS F 2 structure 
functions calculated at Q 2 = 10 GeV 2 . The tar- 
gets are from the top to the bottom: proton, neu- 
tron and their average isoscalar one. The plots 
for the DIS part are done with GRV94 PDFs |7|. 
RS structure functions are yet not re-scaled by 
l7r functions. Similar plots for F 2 at Q 2 — 0.5, 1 
GeV 2 were shown before in [H] . We see that in the 
case of proton and isoscalar target typical man- 
ifestation of local quark-hadron duality: the A 
peak slides along the DIS curve. The duality does 
not seem to apply neither to other resonances nor 
to neutron target even in the case of the promi- 
nent A resonance. 

In Fig. 21 we show the role played by l7r func- 
tions. On the same plot both Rein-Sehgal and l7r 
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Figure 3. Manifestation of quark-hadron duality for 
Fi structure function for CC neutrino reactions on 
proton, neutron and isoscalar target. F2 for RS model 
is calculated at Q 2 = 0.4, 1, 2 GeV 2 and for DIS at 
10 GeV 2 . 
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Figure 4. The comparison of F2 structure func- 
tions for the RS model, re-scaled RS model and the 
DIS. The structure functions are calculated for CC 
^-neutron scattering at Q 2 =l GeV 2 . 



function re-scaled RES F2 structure functions for 
CC reaction on neutron are presented. We see 
that the modifications apply mostly to the region 
of W close to 2 GeV. 

Fig.[5]presents the functions defined in Eq.0for 
three different integration regions i.e. for three 
definitions of the resonance region: W £ (M + 
m w ,W ma x) where W max = 1.6, 1.8, 2 GeV. We 
see that the choice W m ax = 2 GeV makes the 
functions slowly varying in the wide region of Q 2 . 
The behavior for small values of Q 2 < 0.5 GeV 2 
is very different in agreement with the predictions 
made in ||] . For the proton target the best choice 
is W max = 2 GeV while for other two targets it is 
preferable to choice the resonance region as more 
confined at the price of significant variations with 
Q 2 . The choice W max = 2 GeV is a natural one 
for the RS model. But it is suggested in ^IJI thai 
the RS model underestimates the cross section for 
W > 1.7 GeV. We conclude that with the choice 
W max — 2 GeV the duality holds very well for 
the proton and badly for other targets. It seems 
to be difficult to have duality for all the targets 
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simultaneously. 
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Figure 5. Function TZ 2 (Q 2 ;Q 2 DIS = 10 GeV 2 ) 
(see Eq. |3J for CC reaction on proton, neutron and 
isoscalar targets with W m a X = 1-6, 1.8, 2.0 GeV. 



